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Abstract 

We construct a one-loop effective metric describing the evaporation 
phase of a Schwarzschild black hole in a spherically symmetric null- 
dust model. This is achieved by quantising the Vaidya solution and by 
chosing a time dependent quantum state. This state describes a black 
hole which is initially in thermal equilibrium and then the equilibrium 
is switched off, so that the black hole starts to evaporate, shrinking to a 
zero radius in a finite proper time. The naked singularity appears, and 
the Hawking flux diverges at the end-point. However, a static metric 
can be imposed in the future of the end-point. Although this end-state 
metric cannot be determined within our construction, we show that it 
cannot be a flat metric. 
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1 Introduction 



The two-dimensional dilaton gravity models turned out to be very useful 
toy models of black hole formation and evaporation |JlJ]. Their relevance for 
4d black holes comes from the fact that the spherically symmetric scalar 
field collapse can be described by a 2d dilaton gravity action 
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where G is the Newton constant and the 4d line element ds 4 is related to 
the 2d line element ds by 



dsA 



(1.2) 



R is the 2d scalar curvature associated with the 2d metric g^, <f) is the 
dilaton field and fa are N matter scalar fields. These fields depend only 
on time t and radial coordinate r, while the angular dependence resides in 
d 2 . The spherically symmetric collapse was studied by several authors 
and the problem of determining a semiclassical metric which includes the 
back-reaction of the Hawking radiation is still unsolved. This is related to 
the fact that the classical equations of motion are not solvable. In contrast 
to this, a string theory inspired 2d dilaton gravity model []3T] 
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is classicaly solvable, and its solution describes a formation of a 2d black 
hole. The quantization of ( |1.3[ ) is made simpler by the fact that the matter 
fields propagate freely [||, 0, |, g, 0. 0. 01. so th at the one-loop [0, 
15 , [Tot ], and the two-loop [ Jl6l ] effective metrics were obtained. Therefore 



one can study analytically the back-reaction effects in this 2d model. 

In this paper we are going to study a more realistic 2d dilaton gravity 
model, which will have some of the nice features of fll.3Q but it is going to 
describe a 4d black hole. We will study 
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whose 4d interpretation is that of a self-gravitating spherically symmetric 
null-dust cloud. When compared to the action (]1.1]), one notices that a 
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simplifying feature of ( |1.4[ ) is that the matter fields do not couple to the 
dilaton, so that one obtains free-field matter equations of motion in the 



conformal gauge, in analogy to the action (1.3). This means that the task 



of determining the back-reaction in the model ( |1.4| ) is going to be simpler 
than in the model < \l.ty . Still, the quantization of ( |1.4| ) is complicated by 
the fact that the general solution of the equations of motion is not known. 
However, in a special case when 

ds 2 = _( 1 _ 2m(u)\ ^ + ^ (Lg) 



r 

where exp(— </>) = r, G = 1, the classical solution for fa = fi(v) is given by 



dm(v) 1 / dfi\ 2 



dv 2 ~ \ dv 

i=i 



which is the well-known Vaidya solution pTTp . It describes a collapse of 
a spherically symmetric null-dust cloud. The equations ( |1.5| ) and ([L6|) 
will be the starting point for our quantization procedure, from which we 
will determine an effective metric describing the one-loop back-reaction 



effects. Note that the back-reaction effects in the model ( |1.4| ) have been 
studied in [Q, |5p where the one-loop back-reaction has been modeled by 
adding a Polyakov-Liouville term to the action ( |1.4j ). However, the resulting 
equations of motion are not solvable, and only a numerical study has been 
done. 

In this paper we will perform an operator quantization of the equa- 
tions dl.5D and dL6l), so that an explicit expression for a one-loop effective 
metric will be obtained. This metric will describe the evaporation of a 
Schwarzschild black hole which was initially in a thermal equilibrium 
state. This is achieved by using a quantization formalism developed in 
[0, [lltl, and by using the idea of thermal bath removal Q18Q, which was 



developed in the case of the 2d model ( |1.3[ ). We first show that the idea 
of thermal bath removal can be naturally formulated in the operator for- 
malism, where it corresponds to the introduction of a time dependence in 
the Heisenberg quantum state of the system. This time dependence can 
be attributed to the external forces which switch off the thermal equilib- 



rium. Then we implement this idea to the model ( |1.4j ) and obtain a one-loop 
metric whose properties we study. 
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2 Operator formalism and thermal bath removal in 
the CGHS model 



The general solution of the classical equations of motion following from 
( |1.2| ) in the conformal gauge ds 2 = —e 2p dx + dx~ are (up to constant shifts 
in the x ± coordinates) 

e ~2p = e -24> = _x 2 x + x- - F + - F_ + ^ , (2.1) 

A 

fi = f i+ (x + ) + /<_ (x-) , (2.2) 



where 



.± 



F± = J J T±±(x+) , (2.3) 

and T±± = \{d±f) 2 . M/X is an integration constant, and the residual 
conformal invariance has been fixed by the so-called "Kruskal" gauge p = <fi. 



From the expressions ( |2.1] ), ( |2.2| ) and ( |2.3| ) it is clear that the indepen- 
dent degrees of freedom are those of the matter fields. Therefore a reduced 
phase space (RPS) quantization should give a physical Hilbert space which 
coincides with the Hilbert space of massless scalar fields propagating on a 
flat background [[7], [10| ]. As far as the problem of diffeomorphism anomalies 
is concerned, it is formally avoided in the RPS quantization, although it 
may be hidden in the non-covariant form of the gauge-fixed theory. How- 
ever, an anomaly-free Dirac quantization of the CGHS model gives the 



same physical Hilbert space as the RPS quantization [ ]13| ], which guran- 
tees the diffeomorphism invariance of the RPS results. The dynamics is 
generated by the free-field hamiltonian of N massless scalar fields, and 
therefore the quantum evolution is unitary. The effective metric is de- 
termined by ds 2 = — (e 2p ^j dx + dx~ . The effective conformal factor can be 
evaluated perturbatively by using a matter-loop expansion plOj , [Dp, so that 
at the one-loop order one obtains 

e- 2 ^ = (Vol e~ 2 P |Vo) = -\ 2 x + x~ - (F+) - + y , (2.4) 



where now F± are operator valued expressions ( fO| ) in the Heisenberg pic- 
ture. The initial state |V>o) can be chosen to be a coherent state e A+ \0 a +) <g> 
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|0 CT ), corresponding to a left-moving pulse of matter, where a ± are asymp- 
totically flat dilaton vacuum coordinates [\x ± = ±e ±Aa± ) and |0 CT + ) (g> \0 a - ) 
is the corresponding vacuum. If the normal ordering in T± is chosen to be 
with respect to the Kruskal vacuum |0 a 

N 1 
2 
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where we have chosen the conventional normalization of the flux [ ]19[ ], so 
that there is no a factor of 1/n in ( |1.3| ). The expression ( |2.4| ) then gives 
an evaporating black hole solution corresponding to the one-loop effective 



action of Qlq] 

S 



So 
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Note that one can consider a different process of black hole evaporation, 



if a different initial state is chosen. Instead of the boundary conditions ( |2.5| ), 
which correspond to the gravitational collapse situation, one can consider 
an evaporation process where initially one has a black hole in thermal 
equilibrium and at x + = x^ the incoming thermal flux is switched off [[I8p. 
This process can be described by the boundary conditions 
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It is not difficult to see that the boundary conditions ( |2.7[ ) correspond to the 
following state |o) 



e(x 
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A novel feature of ( |2.8p is that the Heisenberg state | ) now depends on 
time, which reflects the nature of the new process where an external force 



has to be used in order to do the switching. It follows from ( ^8| ) that 
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6(x + -x + ) , <T CT - CT -> 



NX 2 



(2.9) 



48^ v u 48tt ' 

where the normal ordering in T a ± a ± is with respect to the dilaton vacuum, 
so that the incoming constant thermal flux has been switched off. 
The one-loop solution for x + > x^ is now given by 
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where = -iV/487rA 2 XQ . The asymptotically flat coordinates a ± at the 
future null-infinity are given by 



a+ = a + , e~ x ~ a = e- Xa ~ - A . (2.11) 



The first relation in ( |2.11[ ) is consistent with the choice ( |2.8[ ), since it implies 
that for x + > there is no incoming flux at the past null-infinity, i.e. 

(T a+a+ ) = . 



The second relation in ( |2.11[ ) implies that the Hawking flux is the same as 
the initial thermal flux, i.e. 

NX 2 

PW->=4^- (2.12) 

It is not difficult to see that, due to the Hawking radiation, the apparent 
horizon shrinks and meets the curvature singularity in a finite proper 
time. The evaporating solution can be continuously matched to a static 
solution on the null line x~ = xj nt . This solution coincides with the remnant 



geometry of p5[ ], which appears in the evaporation process initiated by a 
gravitational collapse. Therefore this 2d example confirms the intuition 
that the basic features of the evaporation process do not depend on the way 
how the black hole was created. 



3 One-loop analytic model for Schwarzschild black 
hole evaporation 



Now we apply the idea of thermal bath removal to the model ( 1.4). The main 



problem which appears when trying to apply the RPS operator formalism 



to the theory ( |L4D is that, in contrast to the theory (|1.3|), we do not know 
the general classical solution for an arbitrary matter energy-momentum 
tensor T M „. However, if we want to describe the evaporation process of a 
black hole which is initially in thermal equilibrium and then the incoming 
thermal flux is switched off, the problem becomes simpler. 



We start from the Vaidya solution ( |1.5| ), and in analogy with the 2d case, 
we take the following state |o) 

|o) = 6(v - v) |0y) ® \0 V ) + 9(v - v ) \0 V ) ® \0u) , (3.1) 
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where V = 4Mexp(t;/4M) and U = -4Mexp(-u/4M) are the Kruskal 
coordinates of the initial Schwarzschild black hole and M is its mass. Con- 
sistency then requires that for v < w 

(V) = , (3.2) 

which is satisfied if T^ v is normal ordered with respect to the Hartle- 
Hawking vacuum |0y) ® |0t/). The incoming and the outgoing flux are 
constant for v < vq, and take the value corresponding to the temperature 

T = (81TM)- 1 

(3.3) 



48^(4M) 2 ' 

where denotes the operator obtained by normal ordering with re- 
spect to the asymptotically flat coordinates (u, v). For v > vq, one obtains 

^=- 48^)2 =-^ > ( T ^)=0 , (3.4) 



where the first equation follows from 



(T w ) = (0 V | : T vv : v |0„) = ( — j 2 (0„| : Tyv -v |0„) 



N fdV\ 



24-k \dv J 

(3.5) 

and Dy{v) = ■^r — 1 (^7") 2 ^ s the Schwartzian derivative. The effective 
metric is then obtained by taking the expectation value of the expression 
(|L5|), so that 

ds 2 ) = - (l + - \J3v9(v) - r s ]\ dv 2 + 2dvdr , (3.6) 



where we have used (3.2), (3.4) and (1.6). r s = 2M is the Schwarzschild 



radius and we have set vq = 0. The effective metric ( |37q ) is of the one-loop 
order since it is only a function of (T^) and it does not depend on {T^ u T pa ) 
or on the higher-order energy-momentum tensors correlation functions. 

For v > the metric ( |3~l3| ) represents an evaporating black hole whose 
mass is linearly decreasing with time. Such a metric was previously studied 



in [20], where it was ad hoc postulated and used to describe the evaporation 



phase of a black hole which was created from a vacuum. Consequently, a flat 
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spacetime was chosen for v < 0, instead of the Schwarzschild spacetime. 
The advantage of our approach is that the operator formalism provides 
metrics which are consistent with the boundary conditions. In this way 
one avoids inconsistencies which may appear due to the ad hoc nature 
of the procedure used in []20p. For example, a flat metric was chosen for 
v > r s /(3 in [ffijp, and since the Hawking radiation is produced, one obtains 
a flat spacetime with non-zero energy-momentum tensor. 

The line element (13.61) can be written in the conformal form (we will 



omit the expectation value) 



ds 2 
where 



9(-v) (l - —) dvdu + 6(v)-£- (l + z- 2(3z 2 ) dvdu , (3.7) 
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and the coordinate u is determined from the equation 



1 1 — (3v/r t 



\Z — Z-\ 



|1 - Z/Zj 



(3.8) 



(3.9) 



The constants z±, A± are given by 



z± 



4/3 



A± ~ w i 1 ± vt+w 



(3.10) 



Note that (3.9) can be also viewed as the equation determining r = r(u, v). 



The requirement that the conformal factor in (3.7) is continuous at v = 
gives 



u 



— = -A- log 



r 

— + z- 
r* 



A+ log 



1 + 



r s z+ 



where 



r + r s log 



1 



u 
2 



(3.11) 



(3.12) 



The relations ( p.llD and ( |3.12D determine the function u = u(u). One can 
now check that the incoming thermal flux has been removed for v > 0, 
since v remains the asymptotically flat coordinate at the past null infinity 

u — > —oo. 

The scalar curvature of the effective metric is given by 



R=-o(r a - Pv6(v)) 



(3.13) 
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so that the curvature singularity is at r = 0. The apparent horizon curve 
is determined by d v r = 0, which gives 

r A H = r s -f3v6(v)- (3.14) 



tah decreases as the black hole evaporates, and the curve ( |3.14| ) intersects 
the r = curve at 

Vint = r s /P , Ui n t = oo , (3.15) 

so that for v > v int a naked singularity appears, see Fig. 1. 



Since the metric (3.6) is a one-loop approximation, it is valid only in the 
region where Rip 2 < 1 (since the loop-expansion of the effective metric is in 
(T n ), which is of the order of (Rlp 2 ) n ). This gives that tab > V2lp, which 
is the expected Planck length cutoff. Note that if one defines a dynamical 
black hole mass Mp as Mp = \tah, then 

d -^ = -P/2 , (3.16) 
dv 

which does not correspond to a thermal evaporation mass equation which 
is given by ^f- oc —M^ 2 . However, because (3 is very small (in physical 
units it is given by j3 = ^(mp/Af) 2 where mp is the Planck mass), the 
difference between the non-thermal evaporation ( |3.16D and a thermal one 
will be noticed only when tah < r c ^ |r s . Therefore the evaporation is 
thermal for tah > r c » lp. 

The Hawking flux Th can be calculated from the expression 

N 

T H = (ol T UFUF | ) = (0u\ T UFUF \0 V ) = ~^D Up (U) , (3.17) 

where up is the asymptotically flat coordinate at the future null infinity 
(v = oo). By using ( |3.9[ ) one can show that as v — > oo 

ds 2 ~ _^£± \ z+ _ Z _\~ A -/ A + ( e -u/r a A +dvd ~ (3 _ 18) 



Hence the asymptotically flat coordinate up is given by 

up = -A + r s ( e -«/^ A + - l) . (3.19) 



By using ( |3.9[ ) and the implicit relation u = u(u) defined by ( |3.11| ) and ( |3.12| ) 
one can work out the Hawking flux from ( 3.17| ) 

T H = %\x + z^-^-'^l + x/z+^h + APix + l + x- 1 ) 
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+ 4/3 2 (x 2 + 2x + 4)--^- 



(3.20) 



where x = r(u, v = 0)/r s = r(u, v = 0)/r s . The behavior of the flux is plotted 
in Fig. 2. 

For early times (u —* — oo or x —* oo) T# is close to the initial thermal 
value 

Ttf « T x~ 4/3 « T exp (4/3 log(2/3) + 4/3 2 u/r s ) , (3.21) 



while for the late times (u — » oo or x — > — z_) it diverges at the end-point 

x = -z_ = 1 - 2/3 + 0(/3 2 ) as 

T H «T |x + z_|- 4/3 (l + 8/3) . (3.22) 

This is an expected behaviour because of the presence of the naked singu- 
larity at the end-point, and it is related to the fact that ^f- ^ at the 



end-point [pOQ. However, as the end-point is approached, the higher-order 
loop corrections become relevant, and the one-loop approximation is ex- 
pected to break down, so that the one-loop divergence could be removed by 
the higher-loop corrections. This actually happens in the CGHS case when 
the two loop corrections are taken into account []16|]. Therefore one can 
expect that the higher-loop corrections will remove the naked singularity. 

4 Conclusions 

Note that our metric is a self-consistent semiclassical solution in the sense 
that its Einstein tensor is proportional to (T M!/ ) by construction. However, 
our metric does not satisfy an additional requirement that the Hawking 
flux is finite at the future null infinity [p0[]. In our case this means that the 
higher-order quantum corrections become important near the end-point. In 



the 2d case, the one-loop metric of [15] satisfies the both criteria; however, 



it is defined only in the weak-coupling region of the space-time, i.e. in the 
region where the higher-order quantum corrections can be neglected. Note 
that in our case the flux stays very close to the thermal classical value 
until very late times u. Therefore one could employ the BPP strategy of 
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removing the naked singularity by imposing a strong-coupling boundary 
at U = Uf, « 0, where U = — 4Mexp(— u/4M), and then in the region 
U > Uf,, V > Vi n t try to impose a static metric such that it coincides with 
( p77| ) at U = U b for V > V int . Also note that U = line is tangential to the 
r = curve at V = Vint, so that one has the same situation as in the BPP 
case. The diference now is that the value of the Hawking flux is infinite at 
U = 0, which is problematic. This is avoided by putting the strong-coupling 
boundary at U = Ub < 0. 
When u = u b , then 

ds = —Cij(v)dvdu , 

where C b (v) = 2 (j^) \u b , so that a trivial solution in the region U > Ub, 

V > Vi n t is 

ds = —dudv , 

where dv = Cb(v)dv. However, this is not a good solution because u^u and 
v ^ v, which means that radiation is present in the flat spacetime region 
U > Ub, V > Vint- This is no surprise, because one expects that the end- 
state geometry cannot be a flat space, but it should be an asymptotically flat 
quantum corrected vacuum geometry, and the corresponding one-loop effec- 
tive action must contain additional counterterms to the Polyakov-Liouville 



counterterm, in analogy to the BPP case [jig]. This quantum vacuum ge- 
ometry cannot be determined within our construction. However, it is clear 
how our construction can be extended. One should find a more general class 
of classical solutions than the Vaidya solutions, and then quantize them ac- 
cording to our approach. These classical solutions could be obtained either 



approximately or by using the global symmetries of the theory [21] 
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Figure I: Kruskal diagram of the one-loop geometry. 
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r=0 



Figure II: One-loop Hawking flux. 



14 



